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We have searched for neutrinoless r lepton decays into £ and V", where £ stands for an electron or 
muon, and for a vector meson {(j>, lj, A'*", or p"), using 543 fb^^ of data collected with the 
Belle detector at the KEKB asymmetric-energy e^e^ collider. No excess of signal events over the 
expected background has been observed, and we set upper limits on the branching fractions in the 
range (5.9 — 18) x 10~* at the 90% confidence level. These upper limits include the first results for 
the lij mode as well as new limits that are significantly more restrictive than our previous results 
for the £(j), IK*°, £K'° and £p^ modes. 

PACS numbers: 11.30.Fs, 13.35.Dx, 14.60.Fg 



INTRODUCTION 

In the Standard Model (SM), lepton-flavor- violating 
(LFV) decays of charged leptons are forbidden; even if 
neutrino mixing is taken into account, they are highly 
suppressed. However, LFV is expected to appear in 
many extensions of the SM. Some such models pre- 
dict branching fractions for t LFV decays in the range 
which can be reached at the present 
B-factories. An observation of LFV would provide unam- 
biguous evidence for new physics beyond the SM. 

A search for LFV t~ decays into neutrinoless final 
states with one charged lepton ^^(e^ or ji^) and a vector 
meson was first performed by the CLEO collaboration in 
the £-(l), i-K*°, £-K*° and £-p° final states in 1998; this 
search set upper limits in the range (2.0 — 7.5) x 10~^ 
[^. Later, Belle obtained upper limits in the range 
(2.0 - 7.7) X 10"^ using 158 ib'^ of data. In this paper, 
we report an improved search for LFV t~ decays [5| into 
a charged lepton and a neutral vector meson (F°), where 
V° includes to in addition to the 0, K*°, K*° and p° [f] 
. The analysis is based on a data sample of 543 fb~^, 
corresponding to 4.99 x 10* r-pairs collected with the 
Belle detector Q at the KEKB asymmetric-energy e+e" 
collider Q taken at the T(4S') resonance and 60 MeV 
below it. These results supersede our previous published 
results 0. 

The Belle detector is a large-solid-anglc magnetic spec- 
trometer that consists of a silicon vertex detector, a 50- 
layer central drift chamber, an array of aerogel threshold 
Cherenkov counters, a barrel-like arrangement of timc- 
of- flight scintillation counters, and an electromagnetic 



[t] While preparing this paper, we became aware that the BaBar 
group had also reported on a search for — > £~u> in a 
preprint [^. 



calorimeter comprised of CsI(Tl) crystals located inside 
a superconducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux-return located outside the 
coil is instrumented to detect mesons and to identify 
muons. The detector is described in detail elsewhere [7[. 
Two inner detector configurations were used. A 2.0 cm 
radius beam-pipe and a 3-layer silicon vertex detector 
were used for the first sample of 158 fb^^, while a 1.5 cm 
radius beam-pipe, a 4-layer silicon detector and a small- 
cell inner drift chamber were used to record the remaining 
385 fb-i [l^. 



EVENT SELECTION 

We search for events, in which one r decays to a 
charged lepton and two charged hadrons (3-prong de- 
cay) while the other r decays into one charged particle 
(1-prong decay) and missing neutral particle(s). We re- 
construct (j) candidates from K~^K~, ui from 7r+7r~7r°, 
K*° from iiT+TT", K*'^ from A'^7r+ and p° from tt+tt^. 

The selection criteria described below are optimized 
from studies of Monte Carlo (MC) simulated events and 
the experimental data distributions. The background 
estimation is based on MC simulations of the reaction 
e~^e^ — *■ T+T^ as well as qq continuum and two-photon 
processes. The r+r" sample corresponding to 1524 fb^^ 
is generated using the KKMC code [uj . The MC sam- 
ples of q q an d two-photon processes are produced using 
EvtGen [l^l and AAFH [13[, respectively, in amounts 
corresponding to the luminosity of the experiment. The 
signal MC events are generated by KKMC assuming a 
phase-space distribution for t decay. The detector re- 
sponse is simulated by a GEANT3 [ij] based program. 

The transverse momentum for each charged track is 
required to be larger than 0.06 GeV/c in the barrel region 
(—0.6235 < cos 6* < 0.8332, where 6 is the polar angle 
relative to the direction opposite to that of the incident 
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6+ beam in the laboratory frame) and 0.1 GeV/c in the 
endcap region (-0.8660 < cos6l < -0.6235 and 0.8332 < 
cos9 < 0.9563). The energies of photon candidates are 
required to be larger than 0.1 GeV in both regions. 

To select the signal topology, we require four charged 
tracks in an event with zero net charge, and the total en- 
ergy of charged tracks and photons in the center-of-mass 
(CM) frame to be less than 11 GeV. We also require the 
missing momentum in the laboratory frame to exceed 
0.6 GcV/c, and to point into the detector acceptance 
(—0.8660 < cos 6 < 0.9563). Here the missing momen- 
tum is defined as the difference between the momentum 
of the initial e+e~ system, and the sum of the observed 
momentum vectors. An event is subdivided into 3-prong 
and 1-prong hemispheres with respect to the thrust axis 
calculated from the momenta of all charged tracks and 
photons in the CM frame. These hemispheres are referred 
to as the signal and tag sides, respectively. We allow at 
most two photons on the tag side to take into account 
initial state radiation. To reduce the qq background, not 
more than one photon on the signal side is allowed for 
the £~(/), £~K*°, i~K*^, l~ modes while not more than 
two photons in addition to tt*^ daughters are permitted for 
the l~u: modes. A charged particle of the type x {x — fi, 
e, K or t:) is identified using the likelihood ratio param- 
eter, P^. This is defined as Px = L^/ {J2x ^^)' where 
is the likelihood for particle type x, determined from the 
responses of the relevant detectors [l3| . For muon candi- 
dates on the signal side we require P^ > 0.95 while their 
momentum should be greater than 1.0 GcV/c. The effi- 
ciency for muon identification is 92% with a 1.2% proba- 
bility to misidcntify a pion as a muon. Electrons on the 
signal side arc required to have P^ > 0.9 and momenta 
greater than 0.5 GeV/c. The efficiency for the electron 
identification is 94% while the probability to misidcntify 
a pion as an electron is 0.1%. 

Candidate </> mesons are selected from K^K^ pairs 
with invariant mass in the range 1.01 GeV/c^ < 
Mk+k- < 1.03GeV/c2 (±4cr). For both kaon daughters 
we require Pk > 0.8. To reduce the background from the 
7 — > e^e~ conversions the cut Pg < 0.1 is applied. 

Candidate uj mesons are reconstructed from Tr+Tr^Tr*^ 
with the invariant mass requirement 0.757 GeV/c^ < 
M^+^-^o < 0.808 GeV/c^ (±3cr). A Tr° candidate is 
selected from 7 pairs with invariant mass in the range 
0.11 GeV/c^ < M-yy < 0.15 GeV/c^. In order to improve 
the oj mass resolution, the tt" mass is constrained to its 
world average value of 134.9766 GeV/c^ for the oj mass 
reconstruction. 

Candidate K*'^ and K*^ mesons are selected from 
K^TT^ pairs with invariant mass in the range 
0.827 GeV/c2 < Mk^ < 0.968 GeV/c^ (±3cr), which 
satisfy the condition Pk > 0.8 for the kaon candidate 
and Pe < 0.1 for both daughters. 

Candidate p° mesons are selected from tt+tt" pairs 
with invariant mass in the range 0.478 GeV/c^ < 



1.074 GeV/c^ (±4a-), requiring that the 
daughter pions have Pk < 0.1, Pe < 0.1 and momenta 
greater than 0.5 GeV/c. In addition, for the e^p" 
mode, we require P^j < 0.5 for daughter pions in order to 
reduce the two-photon background from ee — > eefifj.. 

Figures [1] (a)- (d) show the invariant mass distributions 
of the (p, CO, K*^ and candidates for the — > IJ-~4': 
— » p^Lo, T~ — > p~K*'^ and A' p" modes, re- 

spectively. The estimated background distributions agree 
with the data. The main background contribution for the 
— » mode is due to qq events involving (j) mesons. 
For the t~ £^uj mode the dominant background comes 
from T" Ti^LoVr decay with the pion misidentified as 
a Icpton. The t~ — > ■k~ ■k'^ ■k~ decay is one of the main 
background sources for the (rK*^, £~K*^ and 

£~ p'^ modes. In this background source one pion 
is misidentified as a lepton for all modes while for the 
£^ K*^ and £~K*^ modes one additional pion is 
misidentified as a kaon. 




1 1.01 1.021.031.041.05 



0.8 0.9 

M o{C3ev/c'^) 

TCTCTl 



160F 




0.6 0.7 0.8 0.9 1 



1.1 1.2 1.3 
(Gev/c^) 



0.4 0.6 0.8 1 1.2 1.4 1.6 
M (Gev/c') 



FIG. 1: Invariant mass distributions of (a) <f) K^K~ for 
— > p^(fi, (b) u 7r+7r~7r'"' for —>■ p^LO, (c) — » 
K^Tv^ for — > K**^ and (d) — > tt^tt^ for M p" 

in the region 1.5 GeV/c^ < Mtvo < 1.95 GeV/c^ and -0.5 
GeV< AiJ < 0.5 GeV. The points with error bars are data. 
The open histogram shows the expected t^t~ background 
MG while the filled histogram is the sum of qq and two-photon 
MGs. The regions between the vertical lines are selected. 



To reduce the remaining background from t^t^ and 
qq, the events from the triangular area defined by the 
missing momentum, Pmiss (GeV/c), and missing mass 
squared, m^igs ((GeV/c^)^) are selected for further con- 
sideration. These requirements are summarized in Ta- 
ble U and illustrated in Fig. [2] by the two-dimensional 
plots of pmiss (GeV/c) versus m^jg^ ((GeV/c^)^) for the 
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TABLE I: Selection criteria using Pmiss (GeV/c) and ni^jigg 
((GeV/c^)^) wiiere Pmiss is missing momentum and Tn^^^^ is 
missing mass squared. 
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FIG. 2: Pmiss vs. m^iss plots for signal MC and data for the 
T — > ^p'^ mode. The regions between the vertical lines are 
selected. 

^J,^ p'^ mode. 

For the £~uj {£~K*° and £~K*'') mode, we require 
that the opening angle between the lepton and uj {K*°) 
on the signal side in the CM frame, 0f^ (.0™o), sat- 
isfy cos^CM < 0.88 (cos6'C™o < 0.93). To remove two- 
photon background for the eV'^ modes, we further require 
that the opening angle, a, between the direction of the 
total momentum of charged tracks and 7's on the sig- 
nal side and that on the tag side satisfy the condition 
cos a > —0.999 for the e^cp, cos a > —0.996 for the e^ui, 
e-K*°{K*°) and cos a > -0.990 for the e~ p° mode. 

To identify signal r decays, we reconstruct the £V'^ in- 
variant mass, Miyo, and the energy difference in the CM 
frame, AE^ between the sum of energies on the signal 
side and the beam-energy, E'bcam- Signal events should 
concentrate around M^yo = mr and AE ~ 0, where 
TTir is the nominal r mass. For the £uj modes, we used 
the beam-energy constrained mass, Mbc, instead of the 
invariant mass Minv, where Afbc = \/E'^ca.m ~ iPrY^ 
order to improve the mass resolution, which is smeared 
due to the 7 energy resolution. In calculating the r mo- 
mentum Pt , we replace the magnitude of the tt^ momen- 
tum with the value obtained from the beam energy, the 
energies of charged tracks on the signal side and the tt" 



direction measured by the calorimeter. 

The resolutions in AE and Miyo , evaluated using the 
signal MC, are summarized in Table [III We define the 
signal region in the AE — Algyo plane as a ±3(t ellipse. 
In order to avoid biases in the event selection, we blind 
the signal region until the analysis is finalized. 

TABLE II: Resolutions in Mnyo in MeV/c^ and AE in MeV. 
The superscripts low and high indicate the lower and higher 
sides of the peak, respectively. 
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BACKGROUND ESTIMATION 

After all selections, a few events remain in the region 
-lO'^lTCo < ^^^v" < W^fC and -lOa^^ < AE < 
3cr]^'|?, which we define as a background region. This 
region will be used to estimate the expected background 
and is shown in Fig. [3l Table IIIII lists the numbers of 
events in the background region excluding the signal re- 
gion. The background is efficiently suppressed by the 
event selection. The comparison between the data and 
MC shows reasonable agreement for all modes, except for 



TABLE III: Number of events in the background region ex- 
cluding the signal region for data and MC. 
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FIG. 3: Distributions of l^E - Miyo in the ±10ct box for (a) 
^ M^'/'- (t>) ^ (c) — > fi^uj, (d) —> e^oj, 

(e) r- ^_M"A'*», (f) -> e^A'*", (g)_r- ->_M~A*°, (h) 
T — > e A*'', (i) T ^ /i and (j) r e after all 
selections. Dots are data and filled boxes show the signal 
MC. The elliptical area is the 3a signal region. 



the T~ fJ'~P^ mode. For this mode one of the domi- 
nant background sources is qq, which at low multipUcity 
is poorly described by MC. For all other modes we es- 
timate the number of background events in the signal 
ellipse from the data in the background region using the 
following method. 

Since for most of the modes there are single events only 
in the AE — M/yo plane both in the data and MC, we 
first study the background distribution in a sideband re- 
gion larger than the background region (1.5 GeV/c^ < 
Miyo < 1.95 GeV/c^ and -0.5 GeV< AE < 0.5 GeV) 
and find that the distribution of events in it is approxi- 
mately flat in the background region. Therefore we as- 
sume that this distribution is flat inside the background 
region and estimate the expected number of background 
events in the signal region from the number of data events 
in the background region times the ratio of the areas of 
the signal ellipse and background region. If the number 
of data events in the background region is zero, we assign 
an upper limit of 2.44 events at the 90% confidence level. 
The expected number of background events obtained by 
this method is shown in the third column of Table IIVI 

For the fi~uj mode, where the number of events is 
larger, we estimate the background contribution in the 
signal region using the shape of the background MC dis- 
tribution normalized to the data yield in the sideband 
region of the background region. 

For the mode, the background events in FigO 

come mostly from r~ — > 7r~7r"'"7r~z^ decay and qq when 
one of the pions is misidentified as a muon. To estimate 
the background contribution we select a special event 
sample requiring < 0.1 for muon candidates instead of 
> 0.95. The number of expected background events 
is then calculated from the product of the number of 
events with < 0.1 in the signal region and the muon 
fake rate. 



RESULTS 

After unblinding the signal region single events only re- 
main in some modes, see FiglS] The observed number of 
events in the signal region is consistent with the expected 
background. From the numbers of observed events in the 
signal region and the numbers of expected background 
events, listed in the second and third columns of Ta- 
ble llVi respectively, we evaluate the upper limit on the 
number of signal events at the 90% CL, sgo, with sys- 
tematic uncertainties included in the Feldman-Cousins 
method llf 



using the POLE code [l7|. In the cases 
when we give an upper limit of the expected background 
(t — > e~uj, e~K* and e~/9° modes), the number of back- 
ground events is taken to be zero. This results in con- 
servative upper limits. The main systematic uncertain- 
ties on the detection efficiency come from track recon- 
struction (1.0% per track), electron identification (2.2%), 
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muon identification (2.0%), kaon/pion separation (1.4% 
for (j) reconstruction, 1.1% for K*^ and 1.5% for 7r° 
reconstruction (4.0%), statistics of the signal MC (1.3% 
for 0.7% for 0.6% for i-K*° and £-/f*o, 0.5% 
for and 0.6% for e~ p^) and uncertainties in the 

branching fractions for (f) K^K^ and oj tt+tt^tt^ 
(1.2% and 0.8%). The uncertainty in the number of t- 
pair events mainly comes from the luminosity measure- 
ment (1.6%). 

The upper limits on the branching fractions, B, are cal- 
culated asB < 2^^, where Nrr = 4.99 x 10^, is the total 
number of the r-pairs produced and e is the signal efB- 
ciency including the branching fractions of </) — > K^K^, 
u) 7r+7r-7r°, K*° K+tt" and p° tt+tt" The 
resulting upper limits on the branching fractions are sum- 
marized in Table ITVl 

TABLE IV: Summary of the number of observed events A'obs, 
the number of expected background events A'exp, detection 
efficiency e, total systematic error Ae/e, 90% CL upper limit 
of the number of signal events sgo and 90% CL upper limit of 
the branching fractions B. 



Mode 




A^oxp 


e 


Ae/e 


S90 


UL on B 


r 








(%) 


(%) 




(90% CL) 






1 


0.17 ±0.12 


3.14 


5.2 


4.17 


1.3 X 10"^ 


e~ 


<A 





0.18 ±0.12 


3.10 


5.3 


2.27 


7.3 X 10"* 


M" 







0.19 ±0.20 


2.51 


6.3 


2.22 


8.9 X 10"* 


e~ 


Ul 


1 


< 0.24 


2.46 


6.3 


4.34 


1.8 X 10"'' 




-K*° 





0.26 ±0.15 


3.71 


4.8 


2.20 


5.9 X 10"* 


e~ 


K*° 





0.08 ± 0.08 


3.04 


4.9 


2.35 


7.8 X 10"* 






1 


0.17 ±0.12 


4.02 


4.8 


4.14 


1.0 X 10"'' 


e~ 


K*° 





< 0.17 


3.21 


4.9 


2.45 


7.7 X 10"* 




-p" 


1 


1.04 ± 0.28 


4.89 


4.9 


3.34 


6.8 X 10~* 




p' 





< 0.17 


3.94 


5.1 


2.46 


6.3 X 10"* 



efficiency is improved, e.g. it increases by a factor of 2.8 
for the and 2.5 for the e~(j) mode. The new upper 
limits can be used to constrain the parameter space of 
various scenarios beyond the SM. 
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SUMMARY 

We have searched for the LFV decays r~ — ^ l~oj, 
£-K*°, £-K*° and using 543 ih~^ of data obtained 
in the Belle experiment. No evidence for a signal is ob- 
served, and the upper limits on the branching fractions 
are set in the range (5.9 - 18) x 10"^ at the 90% CL. 
This analysis is the first search for the t" £~uj mode. 
The results for the r" ^ ^"0, £-K*°, £- K*° and £- p° 
modes are 3 — 10 times more restrictive than our previ- 
ous results obtained using 158 fb~^ of data. The sen- 
sitivity improvement comes from a factor of 3.4-times 
larger statistics and an optimized analysis. In particu- 
lar, we have improved the conditions on pmiss and m'^-^igg 
to reduce rr and qq background as well as those on the 
opening angle a to reduce two-photon background. As a 
result, better background suppression is achieved and the 
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